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ABSTRACT 

In previons stndies, it has been shown that the long term time average jet 
power, Q, is correlated with the spectral index in the extreme nltraviolet (EUV), 
o^Euv (dehned by ~ i/-<^buv compnted between 700A and llOOA). Larger 
Q tends to decrease the EUV emission. This is a cnrions relationship becanse 
it connects a long term average over ~ 10® years with an instantaneons mea- 
snrement of the EUV. The EUV appears to be emitted adjacent to the central 
snpermassive black hole and the most straightforward explanation of the correla¬ 
tion is that the EUV emitting region interacts in real time with the jet lannching 
mechanism. Alternatively stated, the Q - oeuv correlation is a manifestation 
of a contemporaneons (real time) jet power, Q{t), correlation with aEuv- In 
order to explore this possibility, this paper considers the time variability of the 
strong radio jet of the qnasar 1442-1-101 that is not aberrated by strong Doppler 
enhancement. This high redshift (z = 3.55) qnasar is nniqnely snited for this 
endeavor as the EUV is redshifted into the optical observing window allowing for 
convenient monitoring. More importantly, it is bright enongh to be seen throngh 
the Lyman forest and its radio flux is strong enough that it has been monitored 
frequently. Quasi-simultaneous monitoring (Eve epochs spanning ~ 40 years) 
show that increases in Q(t) correspond to decreases in the EUV as expected. 
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Introduction 


Powerful relativistic radio jets have been associated with the most intrinsically lumi¬ 
nous ultraviolet emitters in the Universe, namely quasars, since their discovery over 50 years 
ago. Only a small fraction of all quasars (~ 10%) have strong jets (radio loud quasars 
RLQs), and even fewer have radio lobes on sup er-galactic scales, ~ 1.7% of all quasars 
have such extended structure ( deVries et ah 20od) . However, until recently there has been 
no observations of the jet launching region, so there has been much speculation as to the 
mec hanism that creates these jets. The literature is hlled with numerical models and theo - 
ries f Lovelace 1976 : Blandford and Znaiek 1977 : Blandford and Pavne 1982 : Punslv 2008 1. 
There is no discriminant for authenticity due to a lack of direct measurement of the jet 
launching region. However, this has changed recently with the examination of extreme ul¬ 
traviolet (EUV) spectra shortward of the peak of the spectral energy distribution (SED) at 
llOOA. The quasar luminosity is widely believed to arise from the vis cous dissipation of tur¬ 
bulence driven by the differential r otational shearing of accreting gas ( Lvnden-Bell and ReesI 


197ll: IShakura and SnnvaevI Il973h . In numerical and theoretical models, the highest fre¬ 


quency optically thick thermal emission arises from the i nnermost region of the accretion 
flow and its frequency is shortward of the peak of the SED (jZhu et al.ll2012l) . The EUV spec¬ 
trum beyond the peak of the SED of quasars is the putative Wien tail of the emission of the 
i nnermost thermal cq r nponent of th e accretion flow that is adjacent to the central black hole 
flMarshall et al.l 119971: IPunslvl 1201411 . If Q is the long term (~ 10® years) time averaged jet 
power (as determined from radio lobe emission and morphology), and Lhni i s the bolometric 


thermal emission associated with the accretion flow, it was shown in iPnnslvl (12014 l2015b|l 
that jet efficiency, Q/Lboi, (which depends on a long term average) was correlated with the 
dehcit of EUV emission in RLQs relative to radio quiet quasars quantihed by aEuv (fhe 
flux density scales as and is computed between llOOA and 700A). This is the 

fundamental correlation and a partial c orrelation anal ysis shows that the correlation that 
exists between Q and aEuv is spurious (IPunslvl l2015bl) . This provides the hrst connection 
between jet power and an observable from a region that is likely contiguous or coincident 
with the jet launching region in quasars. 

This is a curious circumstance since the plasma in the radio lobes that is used to de¬ 
termine Q was ejected from the central engine ~ 10® years before the EUV emitting gas 
reached the environs of the central black hole. Two parameters, L^oi and oeuv are “real 
time” diagnostics of the quasar at the time of observation and Q is a long term time aver¬ 
age. Why should these parameters be connected in what is likely a time variable system? 
It was concluded based on a statistical argument that the most logical explanation is that 
there i s an underlying real time c onnection between the instantaneous jet power, Q{t), and 
a EUV (iPnnslv and Marzianil 120151 ). It is too much of a coincidence that two highly dynamic 
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elements emanating from a common region are strongly correlated if there is no real time 
connection. The scatter seen in the correlation (see Figure 1) includes the variation of Lhoi 
and aEuv from their long term time average value and the degree of scatter indicates that 
in general these variations are modest. 


In this paper, a search is initiated for the real time connection between Q{t)/Lboi and 
aEuv- This is a daunting task since it is not trivial to hnd an estimator that scales with Q{t) 
that is reliable and there is simply not enough observing time on the Hubble Space Telescope 
(HST) to monitor quasar continua (spectra are needed to excise the emission lines) in the 
EUV. Sufficient monitoring requires observing a source with earth-based telescopes which 
requires 2 ; > 3 in order to get the EUV redshifted into the U or B bands. This is not ideal due 
to the large flux attenuation by the intervening Lya forest. Firstly, the quasar must have an 
incredibly large intrinsic EUV luminosity in order to shine bright enough to be seen easily 
through the absorbing screen. The second requirement is that the radio luminosity should 
be high enough so that the object appears frequently in survey work so that there is an 
ample database to cull through for quasi-simultaneous observations. The third requirement 
is almost mutually exclusive from the hrst two. In order to monitor Q{t) by radio flux, the 
source cannot suffer from Doppler aberration. This is necessary since small variations in 
geometry would appear as signihcant changes in the obse rved flux and this would p rovide 
a false indicator of changes in the intrinsic jet luminosity flLind and Blandfordlll985l) . Most 
of the radio sources at 2 ; > 3 in the older radio catalogs have sufficient flux density because 
they are Doppler enhanced. Thus, we have a very strong set of restrictions due to the high 
redshift, luminosity and lack of Doppler beaming: the source needs an enormous intrinsic 
radio luminosity with no signihcant Doppler beaming (i.e. ~ IJy at z > 3 with no Doppler 
beaming). Do such sources exist? They seem to in the form o f gigahertz p eak radio sources 
(GPS). Some of these at high redshift are incredibly luminous fjO’Dealll998l ). We argue below 
that the quasar 1442-1-101 at z = 3.55 is a suitable candidate. 


2. Source Selection 

The high redshift quasar, 1442-1-101, is a powerful radio source that appeared in many 
of the early radio surveys. The spectrum peaks at ~ 1 GHz at 2.5 Jy. Therefore, it was very 
surprising that the source was at extremely high redshift. It was also surprising that unlike 
other high redshift sources that appear in radio surveys, it is not blazar-like in that it is not 
highly variable at radio frequencies. 








Ana- 


4 


EUV Spectral Index Versus Jet Efficiency 
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Fig. 1.— The correlatio n of aEuv versus the jet efficiency, Q/Lboi fjPunslvl l2015bl: 


Puuslv and Marzianill2015l) . More powerful jets have depressed EUV emission. 




















































FA(units of 10'^^ ergs/s/cm^/A) 


- 5 - 


SDSS DR10 Spectrum of 1442+101 


70 



0 ^^^^^^^^^^^^- 

3600 4100 4600 51 00 5600 6100 6600 71 00 7600 8100 8600 91 00 9600 10100 


Wavelengt\h 


Fig. 2. 


Note the ample flux shortward of Lya in the SDSS DR 10 spectrum of 1442+101. 
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2.1. Doppler Beaming 

It is argued here that the radio emission at frequencies below 10 GHz is not Doppler 
beamed in 1442+101. A list of evidence of Doppler beaming occurs below. For each item in 
the list, the actual circumstance observed in 1442+101 is listed after the colon. 


1. Extreme radio variability: Not evident in the large dataset co r npiled in the NASA Ex - 
tragalactic Database or the monitoring in iTingav et all (l2003l) : iMinealiev et al.l (1201211 . 
This is indicative of non-blazar radio emission. 

2. Superluminal motion on 30 pc - 100 pc scales: There is no measurable change in 
position of the three components seen at 1.6 GHz VLBI imag es in 1989 and i n the 
2.3 GHz VLBI images from 1 999 on scales of 10 pc - 100 pc ( Dallacasa et ah 1995 


Pushkarev and Kovalev 11201211 . No blazar-like pc scale evolution has been observed. 


Flat spectrum radio core to > 30 GHz in the quasar rest frame: The spectrum is 
very steep with a spectral index a > 1 ab ove 30 GHz in the quasar rest frame 

The spectrum continues to steepen 


Kovalev et ah 

1999: 

Mantovani et ah 

2009) 


towards 1000 GHz in the quasar rest frame fISteppe et all 1199511 . Thus, there is no 


“buried” flat spectrum core embedded within the source of the gigahertz peaked emis¬ 
sion. There is insignificant blazar-like synchrotron self-absorption at high frequency. 

4. Large optical variability: Low optical variability is reported ( Pica et ^Il988 : Smith et af 


19931 1. The optical variability is similar to radio quiet quasars and thus not typical of 


blazars. 

5. Large optical polarization: No measurements exist, 
activity. 


Inconclusive evidence of blazar 


The preponderance of evidence indicates that Doppler beaming does not affect the emission 
from 1442+101 signihcantly. 


2.2. Estimator for Q{t) 

If a source is not Doppler beamed then it seems reasonable to choose the optically thin 
radio luminosity created in the most compact regions of the source as a quantity that scales 
with jet power. This emanates inside of any working surface (a region of intense jet dissipa¬ 
tion such as an interaction with a dense cloud) and without synchrotron - self absorption, we 
are looking directly at the synchrotron emission that directly reflects dynamical elements. 
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the energy of electrons and the energy of the magnetic held in the jet. It might not be a 
linear relationship, bnt increases in the optically thin radio emission from the innermost jet 
should be positively correlated with the jet power. The quasar 1442+101 is one of the most 
compact radio sources know n (Ivan Breugel et al.lll984l) . There is no emission, even in deep 
observations, on kpc scales ( Murphy et ah EqQsI) . Based on VLBI (very long baseline inter¬ 
ferometry) images of GPS quasars at lower redshift (higher spatial resolution than is possible 
with 1442+101 observat ions), the radio ernission is hig hly resolved with a small fraction in 
the core unlike blazars (IStanghellini et al.l 119971 1200 ll ). Most of the emission especially at 
lower frequencies comes from knots in the jet > 10 pc from the core, the putative working 
surface. This is the primary source of the optically thick low frequency emission in most 
GPS quasars and likely 1442+101 as well. To estimate Q{t), it is not appropriate to sample 
the optically thick radio emission from the “working surface” that depends on dissipation 
that is determined by dynamics and physical conditions external to the jet. Above 8.4 GHz 
(observer fr ame) the spectral index is a > 1 (optically thin) and this frequenc y is sampled 
frequently ( Kovalev et ah 1999 : Mantovani et ah 2009 : Mingaliev et ah 2012 ). Thus, the 
8.4 GHz flux density should be an indicator of jet strength well inside the working surface. 
VLBI imaging at 8.6 GHz indicate s that the majority of the 8.4 G Hz flux is unresolved inside 
of 50 light years from the source (jPnshkarev and Kovalev II2012I ). This is the surrogate to 
monitor jet power given the insufficient monitoring at higher frequencies. 


2.3. EUV Penetration of the Lyman Forest 


The most startling aspect of 1442+101 is that the EUV emission is so bright that it 
shines through a dense forest of Lya absorbing clouds. The spectrum in Figure 2 indicates 
that there is more than adequate luminosity to sample the spectrum at Aq = 4000 A (observer 
frame) corresponding to Ae = 880 A in the quasar rest frame. The flux measured is not 
the intrinsically emitted flux, but that attenuated by the Lyman forest. About hve or 
six absorbers have been resolved spectroscopically in a 6Xo = 100 A window centered at 
Ao = 4000 A ( Baldwin et aLlll974l: iBarthel et aLlll990l) . The intervening gas is comprised 
of many small absorbing clouds, so one does not expect signihcant time variation in the 
total flux integrated over a 100 A window due to changes in the total absorbing column. 
By m easuring the EUV continuum at Ae = 880 ± 11 Awe avoid contamination by emission 
lines ( Teller et al.l 12002 ). The 100 A “smoothing” window minimizes any affects caused by 
spectral resolution differences between observations and the narrow absorption lines, allowing 
for a consistent method of determining the observed flux density at Ao = 4000 A. 






































3. Results 


There are far more radio observations than calibrated optical spectra down to Ao = 
4000A. The optical and radio variability is small fjPica et al.lll988l: iTingav et a]|l2003[) . Due 
to the paucity of optical data sampling it thus seems reasonable (but not ideal) to consider 
observations within 300 days as quasi-simultaneous. Due to cosmological redshifting, this 
corresponds to about 66 days in the quasar rest frame. To put this in perspective, we estimate 
the size of the putative active EUV/radio emitting region in 1442+101 as follows. Based 
on standard bolometric corrections, Figure 2 and archival IR (rest frame optical) data from 
the NASA Extragalactic Database, we estimate L^oi ~ 4 — 5 x 10"^'^ ergs/sec, ignoring the 
contribution fr o m reprocessed emissio n (avoiding double counting) in the molecular clouds 
fjPnnslvl l2015bl: iDavis and Laorl l201ll) . If the accretion flow is radiating at < 50% of the 
Eddington rate, the black hole mass would be M^h > 8 x IO^Mq. For a rapidly spinning 
black hole, 66 days corresponds to < 2 Keplerian orbital periods at a distance of 5Rg {Rg > 
1.1 X 10^® cm is the black hole radius for a rapidly rotating black hole) from the central 
black hole in the equatorial plane. This also corresponds to < 4 Keplerian orbital periods at 
a distance of 3Rg from the central black hole in the equatorial plane. These are reasonable 
estimates for the location of the EU V emitting region from the variability of the EUV and 


numerical models of accretion flows (iMarshall et al.lll997l: iHaba et al 1120031: IZhu et al.ll2012 


Punslvi I2014 [2015b|). Thus, 300 days is not that unreasonable of a condition for quasi¬ 
simultaneity in the context of jet production and EUV emission. 

Table 1 shows the data for the 5 quasi-simultaneous epochs that are spread out over 
~ 40 years (1973 - 2012). Table 1 lists the date of observations (column 1), the 8.46 GHz 
flux density (column 2), the flux density at Aq = 4000A (column 3), the telescopes (column 
4), the difference in the observation date from the optical and radio in days (column 5) and 
the references (columns 6 and 7). For observations that are at a slightly different frequency 
than 8.46 GHz (~ 100 - 400 MHz), a spectral index of 1 is assume d based on the MJD 


47266 VLA observations at 8.21 GHz and 8.66 GHz and archival data (Kovalev et al. 


Minealiev et al. 

2012; 

Mantovani et al. 

2009 

). We note that the observation in 

Barthel et al. 


1999 


of all the epochs. Issues with the flux calibration have been noted based on the atmospheric 
dispersion correction for the iBarthel et al.l fll990l) sample flGorbiru Il99ll) . For SDSS DRIO 


quasar data, such as the MJD 55976 observation in Table 1, there is a calibration issue in the 
blue part of the spectrum. Specihcally, in order to increase spectral sensitivity for quasars, 
the hole drilled over the aperture is centered on the displacement associated with 4000A. 
However, the calibration of the star has a hole centered about the displacement associated 
with 5400A. This inconsistency artihcially raises the flux level in the blue that we estimate 
as a 7.2% excess for quasar light at 4000A (from the airmass of 1.1 and seeing of 1.4 arcsec). 






















































9 


The Anticorrelation of the EUV Flux and the Radio Flux 
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Fig. 3.— The top frame shows the light curves for EUV flux density and 8.46 GHz flux 
density sampled at quasi-simultaneous observation epochs. The bottom frame shows the 
anti-correlation of the EUV flux density and the 8.46 GHz flux density. The squares indicate 
the observational epoch that is not in the “gold sample” (not a VLA radio observation) 
described in the text. 
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Table 1: Quasi-simultaneous EUV and Radio Variability of 1442+101 


Observation 
Date (MJD) 
Radio/ Optical 

Flux Density 
8.46 GHz (mJy) 

Flux Density 

Ao = 4000 A(Ae = 880 A) 
(10-17 

(ergs / sec / cm^ / A) 

Telescope 

Radio/ Optical 

Delta®^ 

(Days) 

42126/ 42189 

598 ± 42 

23.4 ± 2.4^ 

Parkes 64m/ Lick 3m 

297 

47266/ 47320 

709 ± 18'^ 

15.3 + 2.3® 

VLA/ Hale 5m 

54 

48415/ 48427 

704 ± 18'^ 

16.4 + 1.2 1 

VLA/ Lick 3m 

12 

53940/ 53827 

663 ± 16 ^ 

21.3 + 1.3*1 

VLA/ SDSS DR7 

113 

55834/ 55976 

659 ± 16 ^ 

20.5 + 1.4^ 

VLA/ SDSS DRIO 

142 


“Number of days between radio and optical observations 


^ Shimmins and Bolton (1981) 


' Baldwin et al 

(19741 

‘^This paper 

Barthel et al. 

(199nll 

t Lyons et al. ( 

1995 

) 


Piecewise Continuum Spectra Comparison 



Rest Frame Wavelength (Angstroms) 


Fig. 4.— The piecewise continuum estimates in the far UV and EUV of the different epochs 
labeled by their 8.46 GHz flux density. The squares indicate the observational epoch that is 
not in the “gold sample” (not a VLA radio observation) described in the text. The continuum 
excess of weak jet states relative to strong jet states is largest in the EUV. 
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If there is a similar underestimate of the standard star flux at 4000A, the total estimated 
increase in the DRIO quasar spectral flux is 14.4% at 4000A. As a check of this calculation, 
we note that this agrees with the offset seen between the photometric g magnitude and 
the synthetic magnitude computed from the spectrum (18.61 and 18.39, respectively). The 
corrected spectrum is plotted in Figure 2. 


Table 1 can be used to asses s the re a l time connection between Q(t)/Lhoi and aEuv- 
Notice that, per the discussions by Punslv f 2015bll . L^oi ~ ^■8Fx^{Xe = llOOA) in the quasar 

rest frame and aEuv = log -^^e(^e-880A) ^ Rooa^ 0 Thus, aEuv varies monotonically with 

FAg(Ae=1100A) 880A 

F\X^e = 880A)/FAAAe = llOOA). If there is a relationship between Q(t)/Lboi and ueuv 
then there is a relationship between Q{t)/Lboi and Fx^{Xe = 880A)/FAe(-^e = llOOA). The 
advantage of the latter relationship is that both quantities have Fx^{Xe = llOOA) in the 
denominator. Thus if one multiplies both quantities by Fx^{Xe = llOOA), one only needs to 
look at the light curves of Q{t) (as traced by the 8.46 GHz flux density) and Fx^{Xo = 4000A) 
(equivalently Fx^{Xe = 880A)) in order to detect a correlation between Q{t)/Lboi and aEuv- 
The motivation for this re-normalization is to remove the undesirable statistical scatter 
caused by dividing both the quantities by Fx^{Xo = SOOOA) (corresponding to Ag = llOOA) 
that adds to each measured quantity considerable uncertainty, but does not add physical 
content. Such an expedience would not be justified in an analysis that comprised multiple 
objects since each object is normalized differently and the normalization is time variable. 


The top frame of Figure 3 is a plot of both 8.46 GHz flux density and Fx^{Xo = 4000A) 
as a function of time. The bottom frame of Figure 3 shows that Q(t) decreases as Fa„(Ao = 
4000A) increases as expected. Even though the number of observations is too small for rig¬ 
orous statistical analysis, we note that Pearson correlation coefficient is -0.949 corresponding 
to a one-sided 99.3% statistical signihcance. If one restricts the analysis to a “gold sample” 
of epochs all observed with the same telescope (the VLA) and within 150 days of the opti¬ 
cal observation (< 2 Keplerian orbits 3Rg from the black hole), one still sees a correlation. 
The two stronger radio states have the weakest Fa„(Ao = 4000A) and vice versa for the two 
weaker radio states. 


In the analysis of HST composite spectra in iPunslvl (1201411 : iTelfer et al.l (120021) . a deficit 
of EUV continuum emission was associated with the radio loudness in quasars. There was no 
dehcit of emission seen in the far UV continuum for RLQs. However, the continua plotted 
in Figure 4 indicate a far UV continuum suppression as the jet gets stronger. The plot 


^Our spectra d o not extend all the way to Ae 
for other quasars ( Punslv 2015bll . Using Ae = 
expedience to see changes in spectral slope. 


= 700 A. A power law fits the continuum well in this range 
880 A to compute the power law should be an adequate 
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shows piecewise continuum spectra estimates in the far UV and EUV at the different epochs 
labeled by their 8.46 GHz flux density. The rest frame points were sampled at wavelengths 
in which line emission seemed minimal in all the spectra. Rest frame wavelengths were 
indicated for ease of comparison to well known emission lines. The observed wavelengths 
are 4000 A, 5000 A, 6200 A, 6700 A, and 7500 A. There are three radio “weaker” jet 
states at 598 mJy, 659 mJy, 663 mJy and two “stronger” jet states at 704 mJy and 709 mJy. 
The UV fluxes are highly clustered at Ae = 1650A, four of the flux density measurements 
are not statistically significantly different. The flux density of the continua associated with 
the weak jet states exceeds that of the continua associated with the strong jet states by an 
ever increasing differential as the the wavelength decreases. Furthermore, the EUV spectra 
from llOOA to 880A are softer (larger aEuv) for the stronger jet states than for the weaker 
jet states as expected from a Q{t)/Lhoi and ueuv correlation. The continuum suppression 
associated with jet power might taper off in the far UV as opposed to an abrupt turn off at 
Ae = llOOA. 


Discussion 


The physical interpretation of the correlation between jet power and the decrement in 
EUV luminos i ty that is illustrated in Figure 1, 3 and 4 was elucidated in two previous studies 
fjPunslv 2014 . 2015b '). In this section, we review this analysis and discussion in order to put 
the results into perspective. The most straightforward explanation of the correlation involves 
jets from magnetic flux in the inner accretion disk where the EUV emission also originates. 
As more poloidal vertical magnetic flux is stored in the innermost accretion disk, the volume 
available for the thermal emitting gas is displaced. More poloidal magnetic flux equates to a 
stronger jet and less optically thick thermal gas equates to a weaker EUV luminosity. Th is 
dynamical conhguration is known as magnetically arrested accretion fjigumenshchev 2008 '). 


Many different dynamical systems have been classihed as magnetically arrested in the 
numerical simulation literature. The relevant numerical simulations for RLQs are the ones in 
which the vertical magnetic flux is distributed in magnetic islands of low plasma density that 
are concentrated in the inner accretion flow. The configuration is far from time stationary. 
The islands are dynamic and buoyantly move through the ram pressure imposed by the dense 
thermal EUV emitti ng plasma by means of a series of Kruskal-Schwarzschild instabilities 
fllgnmenshchevl l2008h . The ram pressure of the accretion flow concentrates magnetic flux 


near the black hole. 
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4.1. A Review of Evidence for Magnetically Arrested Accretion in RLQs 

Ev i dence o f the detailed predictions of magnetically arrested accretion was fonnd in 


Pnnslvl fl2015bl) by exploring a basic model. The hrst prediction is that dne to the ram 
pressnre interaction of the accreting gas with the magnetic islands, the correlation of Q/Li,oi 
and aEuv shonld be stronger than the correlation between Q and aEuv- This was verified 
by a partial correlation analysis that indicated that the correlation of Q/Lboi with oeuv is 
statistically significant and the correlation of Q with oeuv is spnrions. 

The second piece of evidence snpporting magnetically arrested accretion that was fonnd 


m)ie 

m iPnnslvl (l2015b|) is a verification of the specific relationship between the jet power and the 
EUV snppression that arises from the basic model of magnetically arrested accretion. The 
derivation will not be repeated here, bnt the introdnction of some notation is reqnired to 
explain the exact resnlt that was demonstrated. If the model is representative of the physical 
circnmstance depicted in Fignre 1 then the filling fraction of the inner accretion disk with 
magnetic islands, /, is related to the deficit of EUV Inminosity, where the EUV Inminosity 
deficit is evalnated relative to Inminosity of the EUV continnnm of radio qniet qnasars 
(RQQs). In order to compa re the EUV de ficit from object to object, the EUV spectral 
Inminosity was normalized in iPnnslvl fl2015b[) to the spectral Inminosity at the approximate 
peak of the spectral energy distribntion (SED), 


Normalized EUV Spectral Lnminosity = ^ 


L^(A = llOOA) 


( 1 ) 


where Ly{\ = 700A) is the EUV continnnm spectral Inminosity and = llOOA) is the 
spectral Inminosity at the approximate peak of the SED. Defining the EUV deficit of a RLQ 
relative to RQQs reqnire s defining a fid ncial normalized EUV spectral Inminosity for RQQs 
from composite spectra (lPnn.s]vll2015bl ) 


Fidncial Normalized EUV Spectral Lnminosity = ^ 


= llOOA) 


RQQ 


( 2 ) 


where Ly{\ = 700A) is the EUV continnnm spectral Inminosity of the RQQ composite 
spectrnm and T,y(A = llOOA) is the spectral Inminosity at the approximate peak of the SED 
in the RQQ composite spectrnm. The normalized EUV deficit is approximately eqnal to the 
magnetic finx fill fraction, /, of the inner disk (eqnivalently the fractional volnme of EUV 
emitting gas that is displaced by magnetic finx tnbes) in each individnal RLQ 


EUV Deficit = 


1 

■ L^(A = 700A) ■ 

'L^{\ = iiooA)' 



= iiooA). 

. L^(A = 700A) . 

RQQ 


/ 


(3) 
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Since the jet power scales with the square of the poloidal magnetic flux contained within the 
jet base, the s caling t h at is p redicted by the basic model of magnetically arrested accretion 
was shown by IPnnslvl fj2015bl) to be of the form 


1 

■ L^(A = 700A) ■ 

= iiooA)' 

1 2 


_L^{X = llOOA). 

. L^(A = 700A) . 

RQQ 


g/Lboi 


where A is a constant. The best £t to the RLQ data was found to be 


( 4 ) 


g/Tboi = A 


1 

■ L^(A = 700A) ■ 

'Uix = iiooA)' 



.L^(A = llOOA). 

. L^(A = 700A) . 

RQQ 


-I 2.05±0.17 


( 5 ) 


The uncertainty in the exponent includes uncertainty in the hducial RQQ normalized EUV 
spectral luminosity in Equation (2) as well as the statistical scatter in the data. The agree¬ 
ment of the exponent in equations (4) and (5) within these uncertainties provides strong 
evidence that the mechanism that is creating the jet and suppressing the EUV emission is in 
fact similar to magnetically arrested accretion in the innermost accretion disk. It is further 
estimated from the same analysis that / is a few percent for th e RLQs with weaker jets ari d 
~ 50% for the RLQs with the most powerful jets (jPunslvll2ni5bl : IPunsly and Mar7.ianill2m 5h . 


4.2. Review of the EUV Deficit in the Context of Numerical Simulations 


A major shortcoming of 3-D numerical models of accretion onto black holes is that the 
magnetic held topology and poloidal hux distribution is determined by physics beyond the 
assumed single huid ideal magnetohydrodynamic (MHD) approximation. This is critical 
because it has been shown that even small changes in the numerics near the black hole 
will drastically alter the poloidal magnetic hux topology and the entire dynamical scheme 


fjPunslvi l2015aj) . Most prominent amongst these missing features is the microphysics that 


determines the dihusion rate of plasma onto and oh of magnetic held lines and the magnetic 
held reconnection rate. Neither of these are well known in the exotic environment near 
black holes. These physical processes are critical not only for the formation of the magnetic 
i slands, but the time evolution of t he magnetic islands is determined primarily by dihusion 


flTgnmenshchevll2008l: IPimslvll2015al) . Even worse, these dynamical elements occur only as a 


consequence of numerical di husion in moder n simulations in over simplihed ideal MHD single 


huid models of the physics (jPunslvl I2ni5all . Future development is difficult since it is not 


even clear theoretically what the basic principles required for an accurate physical depiction 
would be. Many issues t hat are related to these topics are active areas of invest i gation in 


solar and fusion physics (iRanman et all l2013l: iMalakitI l2009l: iThrelfall et al.l l2012l: lYamada 
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20071). As such it is imperative that observational results such as those presented here are 
needed to guide the course of numerical and theoretical work. 

The observational evidence provided by this study can be used to cull through possible 
numerical and physical scenarios in order to see which ones are viable candidates to represent 
the physical system in a RLQ. The following restricti ons for numer ical work are model 
independent and derive directly from the observations (jPnnslyll2ni5bl ). These can be used 
to guide us toward numerical models that replicate the physics of jet launching in RLQs. 
Empirically, it is determined that both the large scale poloidal magnetic flux at the base 
of the jet and a mechanism that suppresses (but does not eliminate) the EUV emission 
from the innermost accretion flow coexist at the heart of the central engine of RLQs. This 
is too large of a coincidence, some of this magnetic flux must be the same element that 
disrupts, but does not eliminate the innermost accretion flow. This implies that there is 
significant vertical poloidal magnetic flux in the inner accretion flow of RLQs. If a numerical 
effort cannot reproduce this circumstance then perhaps the numerical approximation to 
the relevant physical processes requires further development. Observations are now mature 
enough to lead the numerical work. 

Summarizing the status of the numerical work, it is noted that the magnetically arrested 
accretion discussed above in the intro duction to this sectio n is qualitatively s imilar to the 
3-D numerical simulations explored in iTgnmenshchevI (120081) and IPnnslv et al.l (120091) . They 
produce a £11 fraction, /, in the innermost accretion flow that is consistent with the range 
described below Equation (5). They also provide most of the gross features of time evolution 
required to sup port the jets and suppress the EUV emission in the inn e rmost accretion flow 


(Punslv 


2015blj . By contrast, the simulations in lMcKinnev et al.l (l2012f)dTchekhovskov et ah 


(120111.120121 1 that are heavily seeded with large scale magnetic flux are devoid of magnetic 
islands close to the event horizon. This is evidenced by the claim in iMcKinnev et al.l (120121 1 
that no significant Poynting flux emerges from this region as well as the linked online videos 
of the simulations. The videos show the innermost significant, modest, ma gnetic island con¬ 
centra t ions are located at r > lORn and they are extremely transient. In iMcKinnev et ah 
(120121) : iTchekhovskov et ahl (j201lL l2012l) the fill fraction, / 0 in the innermost accretion 

flow. This is basically a zone of avoidance for vertical poloidal magnetic flux in these simu¬ 
lations. Instead of penetrating the thickness of the equatorial accretion flow, vertically, the 
poloidal field spreads out above and below the accretion disk in radial fans. Thus, these 
simulations are not consistent with the observations of RLQs. Th ese radial fans cq n ipress 
or “choke” the flow. The “magnetically choked accretion flows” of IMcKinnev et al.l (120121) 
heat the innermost accretion flow compressively as part of the magnetic choking process and 
likely increase radiati on from this regio n. This is verified by the recent results from this 
family of simulations (lAvara et al.l 120151) . The inner accretion flow luminosity is claimed to 
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increase significantly over a standard optically thick, geometrically thin accretion disk, the 
opposite of what is observed in RLQs. 


4.3. The EUV Deficit and Jet Formation Across all Scales 


It is tempting to put these results into a broader context. It is possible that all black 
hole related jet phenomena are a consequence of vertical magnetic flux impeding the inner 
accretion flow. In particular, for Galactic black hole binaries (GBH) jet type emission oc¬ 
curs only when the soft X-ray emission (the putative thermal emission from the a ccretion 
disk) is suppressed f R.ussell et al 2011 : Fender et al 2004; Klein-Wolt et al 2002). How- 
ev er, the “coronal” X-ray power law emission is not suppressed and can be very strong 


(jPunslv and RodriguezI 120131: iFuchs et all 1200311 . 


There have also been claims that dips in the coronal X-ray flux (the X -ray power law) 


are coordinated with superluminal ejections in the Seyfert galaxy, 3G120 flMarscher et al. 


2 OO 2 I: iGhatteriee et al.l 120091: iLohhnk et al.l 120131) . If one were to make an analogy between 
the GBHs and 3G120, the superluminal ejections in 3G 120 would have to be the equivalent 
of the discrete superluminal ejections in GBHs. However, superluminal discrete ejections 
in GBHs are characterized by very strong coronal emission in the last hours just before 
ejec tion that reach near historic maxima that can exceed 50% of the Eddington luminos¬ 
ity (iPiinslv and RodriguezI 120131) . During the ejection, the X-ray lum inosity becomes highly 
varia ble with an average value similar to that just before the ejection fjPunslv and Rodriguez 
20131 ). Thus, the claim of coronal X-ray power law luminosity dips in 3G 120 during super¬ 
luminal ejections refers to a different phenomenon. 

The corona X-ray power law emission during superluminal ejections in 3G 120 likely 
emerges from close to the black hole. However, there is no direct evidence that the emission 
from the innermost optically thick accretion flow is suppressed during superluminal ejections. 
The EUV emission from the high frequency tail of the optically thick thermal emission is 
not monitored. As noted above, in GBHs, the X-ray power law luminosity is not correlated 
with the optically thick thermal luminosity. Thus, there is no precedent that justihes the 
assumption that a dip in the X-ray power law equates to a dip in the optically thick thermal 
emission from the inner accretion disk in 3G 120. Thus, the physical state of the inner 
optically thick accretion disk during superluminal ejections is unknown. There is a claim 
of indirect evidence of inner disk movement outward during superlur ninal ejections base d 
on Fe Ka line width arguments derived from three X-ray observations flLohhnk et al.l 120131) . 
However, the interpretation of the X-ray data in t he “h igh state” observed with XMM and 
the inner disk location claimed by iLohhnk et al.l (120131) disagrees with the seminal studies 
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of that data f Ballantvne et al.l 120041 lOele et al.l 12005 ). The original studies of the XMM 
“high state” indicated that the Fe Ka line was produced in cold plasma at > IbRg from 
the central black hole. This is very far from the innermost stable circular orbit where the 
cold plasma should reside if there were evidence of “a complete disk extending down to the 
innermost sta b le circ ular orbit (ISCO) during the XMM-Newton observation” as claimed by 
Lohfink et ah (2013). 


The claim of evidence of an interaction of the jet with the corona during superlumi¬ 
nal ejections in 3C 120 is not well understood, primarily because the corona is not a well 
understood region. There is strong evidence that the corona might be an out-flowing wind 
in Seyfert galaxies flLiu et al.ll2014J: iKharb et al II2015I) . Does the corona envelope the disk 
like a stellar corona, or is it a separate optically thin ra diatively ineffic i ent d isk inside the 


accretion disk, the so-called truncated disk scenario (e.g.. iLohfink et al.l ( 20131) 1? The phys¬ 


ical interpretation of the depressed coronal emission is highly dependent on the model of the 


corona. 


In summary, both GBHs and RLQs show a suppression of the optically thick thermal 
emission when jetted phenomenon occur. The Seyfert galaxy, 3C 120, shows something 
similar, but it is the coronal emission that appears suppressed which does not occur in 
GBHs. It would be a major advance if these three phenomena could be connected. There 
are many poorly understood details that would need to be studied and clarihed further before 
we are at that point. 


5. Conclusion and Future Prospects 

This paper studies archival data on the distant quasar, 1442-1-101. Evidence of a real 
time connection between the jet power and the suppression of EUV emission is found. The 
result is based on only 5 epochs and simultaneity was established only within a 300 day 
interval. These findings motivate the need for a long term monitoring program of this source 
and other high redshift GPS. The radio observations would be improved with multi-frequency 
observations with VLA or ATGA. Knowing the flux density at more than one frequency 
(perhaps 8.4 GHz, 15 GHz and 22 GHz) would allow one to compute the luminosity of the 
optically thin emission and would be less prone to errors induced from slight changes in 
spectral slope that are problematic with a single frequency proxy. Goordinated sampling 
two times a year would be sufficient. Ideally, the source should be monitored with the same 
optical and radio instrument to minimize uncertainties in this difficult measurement. The 
VLA “gold sample” implemented here improves the situation, but there are still multiple 
optical telescopes involved. 
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